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Abstract 

We investigate a description of shape-mixing and shape-transitions using 
collective coordinates. To that end we apply a theory of adiabatic large- 
amplitude motion to a simplified nuclear shell-model, where the approximate 
results can be contrasted with exact diagonalisations. We find excellent agree- 
ment for different regimes, and contrast the results with those from a more 
standard calculation using a quadrupole constraint. We show that the method 
employed in this work selects diabatic (crossing) potential energy curves where 
these are appropriate, and discuss the implications for a microscopic study of 
shape coexistence. 
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I. INTRODUCTION 



In order to describe processes in nuclei involving large excursions from equilibrium, such 
as shape-coexistence or fission, we cannot use the small-amplitude harmonic expansion about 
a stationary mean-field state provided by the random-phase approximation (RPA). Various 
methods exist to deal with large-amplitude motion, many of which are reviewed in Ref . • 
This reference also sets out the basic formalism on which our work is based. 

We have recently embarked on a study into the properties of collective motion in systems 
with pairing. In a first application we have analysed properties of collective motion in a 
semi-microscopic model of nucleons interacting through a pairing force, coupled to a single 
harmonic variable giving a macroscopic description of the remaining degrees of freedom || . 
We have used the local-harmonic approximation, which is equivalent to the generalised valley 
approximation for a single collective coordinate [1] , to analyse the collective motion. To our 
surprise it has turned out that the system automatically selects either diabatic or adiabatic 
collective surfaces according to the strength of the pairing interaction. However, since this 
model is not fully microscopic, we feel that it would be beneficial to study a fully microscopic 
Hamiltonian. This does not mean that we wish to fully ignore the success of the unified 
model by Bohr and Mottelson, which indicates that the semi-macroscopic approach describes 
many nuclear phenomena quite well, but rather that we wish to understand such behaviour 
from a microscopic viewpoint. We thus feel that it is desirable to test the methodology for 
a fully microscopic Hamiltonian which is able to describe nuclear systems from vibrational 
to deformed. 

To this end we investigate the collective motion in a microscopic model which describes 
a system of nucleons interacting through a simplified version of the pairing-plus-quadruple 
force ||. Although the Hamiltonian has a very simple form, we shall see that the model can 
reproduce the qualitative features of many kinds of interesting situations observed in real 
nuclei. In this work we shall concentrate on the spherical-to-deformed transition, and nuclei 
with shape coexistence, where more than one equilibrium shape plays a role. 

In the case of a single- j shell the model Hamiltonian is built from the generators of an 
o(4) algebra, which makes exact diagonalisation feasible. The model has been originally 
developed to describe K 11 = + excitations in deformed nuclei |J, and has also been used 
as a test-bed for various methods used in the calculation of collective excitations such as 
the boson expansion method [|J, the selfconsistent collective coordinate method ||, and a 
semiclassical method || . The model can be generalised to multiple shells, where it has been 
used to investigate shape-coexistence phenomena [J/] . Finally, a similar model has been used 
to study the collective mass parameter in finite superconducting systems ||. 

Although the low-lying spectra in nuclei are mostly dominated by the quadrupole phonon 
(J 71 = 2 + ) excitations, the anharmonicity is very important for many nuclei, especially in a 
shape-transition region, where the nature of the ground-state changes rapidly with particle 
number. For instance, the even-even Sm isotopes show a typical example of the spherical 
to deformed shape transition in which the spectrum shows a strong anharmonicity between 
the spherical (N < 84) and deformed (N > 90) nuclei, especially for 148 ' 150 Sm. These 
phenomena are primarily related to the competition between the monopole and quadrupole 
interactions among the valence particles outside a closed core. The pairing-plus-quadrupole 
model, originally proposed by Bohr and Mottelson, was designed to describe this competition 
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and is quite well able to reproduce the most important aspects of the experimental data (see 
|J and references therein). Later the boson expansion method has been applied to the same 
model (with an additional quadrupole pairing interaction) for the description of the shape 
transition in the Sm isotopes which shows excellent agreement with the experimental 
data. Since the 0(4) model is very similar to the pairing-plus-quadrupole model, it would 
be of significant interest to see whether our method of large amplitude collective motion is 
able to properly describe the shape transition phenomena in this exactly solvable model. 

The importance of shape-coexistence in nuclear physics can be seen from the multitude 
of theoretical approaches and the amount of experimental data as compiled in the review 
paper Ref . |L1[] . An important example can be found in even semimagic Sn and Pb isotopes, 
where the ground states are spherical. However, deformed excited J 71 " = + states have 
been observed at low-excitation energies in many of these nuclei. These excited states are 
regarded as states associated with proton two-particle-two-hole (2p-2h) excitations across 
the closed shell. Using the Nilsson picture, which shows down-sloping single-particle levels 
above the proton closed shell, and up-sloping levels below it, it is possible to assign a con- 
figuration of two particles lying on down-sloping levels and two holes on up-sloping levels. 
The configuration-constrained Nilsson-Strutinsky calculations as performed by Bengtsson 
and Nazarewicz have suggested that the diabatic potential-energy curve obtained by 
switching off the interaction between the 2p-2h and the ground-state (Op-Oh) configuration 
gives more accurate picture than the conventional adiabatic potential energy. This question, 
whether the nuclear collective potential is adiabatic or diabatic, is quite old, and was origi- 
nally raised by Hill and Wheeler ||13|| . We have shown in our previous work p| that, using a 
method to self-consistently determine collective coordinates, the system itself selects either 
an adiabatic or a diabatic collective path according to the properties of the interaction. 
It is our aim to investigate in the 0(4) model whether the method is able to provide us 
with useful information about shape mixing, and to test whether it makes useful predictions 
whether the collective potential energy is diabatic or adiabatic. 

In Sec. |T[ the theory of adiabatic large amplitude collective motion is briefly reviewed. 
A problem peculiar to the model under consideration, related to the special character of the 
zero modes, is dealt with by a prescription to remove the zero-mode degrees of freedom. In 
Sec. |TJ, the 0(4) model in a single- j shell is described, and we discuss the applications of 
our methodology to the model. In Sec. [TV], a generalisation of the 0(4) model to multiple 
shells is introduced. We then investigate the large amplitude collective motion for a set of 
parameters which describe transitional (from vibrational to deformed) nuclei. In Sec. |V], 
we use a proton-neutron form of the Hamiltonian, which can be used to describe shape- 
coexistence. We show how our method selects a diabatic or adiabatic potential energy 
curve. Finally, we give some conclusion and present an outlook in Sec. [VI|. 



II. FORMALISM 



A. Review of methodology for the local harmonic approximation 

We briefly review the local harmonic approximation (LHA), as a method in adiabatic 
large amplitude collective motion (ALACM). A full discussion of the method can be found 
in Ref. [p]] . As usual in our work, we use the convention that the repeated appearance of the 
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same symbol (a, ft, • • ■ ; i, j, • • •) as an upper and lower index denotes a sum over the relevant 
index over all allowed values. We also use the convention that a comma in a lower index 
denotes the derivative with respect to a coordinate, thus F >a = dF/d^ a . 

Our approach to large amplitude collective motion is applicable for classical Hamiltonian 
systems which have a kinetic energy quadratic in momentum only Since most systems, 
especially the mean- field problems of nuclear physics, do not satisfy this requirement, we 
are forced to truncate the full Hamiltonian 7i(£, ix) to second order, 

H ad (Z,7r) = ^B a ^ a np + V(0 , a, ft = 1, ■ ■ ■ , n . (2.1) 

Here the mass tensor B a @ depends in general on the coordinates £ Q and is defined by ex- 
pansion of the Hamiltonian to second order in momenta, 

dn 0,071/3 

Thus all terms of more than quadratic order in momenta are neglected. This is only possible 
when the higher-order terms are small, which is definitely true in the small (zero) velocity 
limit. It is in this sense that the theory may be regarded as an adiabatic theory. The tensor 
B a p, which is defined as the inverse of B a/3 [B ai B 1 p = 5^), plays the role of metric tensor 
in the Riemannian formulation of the local harmonic approximation (LHA) below. 

Collective coordinates x l and intrinsic (non-collective) coordinates x a which are decou- 
pled from each other, are assumed to be reached by making a point transformation, con- 
serving the quadratic nature of Eq. (|2~1~1) , 

z* = /*(£) (i=l,... | J0, (2.3) 
x° = / B (0 (a = K + l,...,n) . (2.4) 

In this section, we use symbols (a, ft, • • •) as indices of the original coordinates, (/i, v, • • •) 
for the new coordinates after the transformation, for collective coordinates and 

(a, &,•••) for intrinsic coordinates. The requirement that the motion is exactly restricted to 
a collective subspace S (defined by x a = 0), yields three conditions, 

B ai = , V a = , (2.5) 
B% = . (2.6) 

These three conditions are only satisfied for exactly decoupled collective motion, a rare 
occurrence indeed. It is therefore practical to combine these three decoupling conditions 
into two sets of conditions, that are satisfied even when no exact collective subspace exists. 
The conditions chosen are those that determine the bottom of a valley in the potential 
landscape, which is, under certain conditions, an approximation to a decoupled subspace. 
The quality of this decoupling can be measured, see below. In this work we prefer to work 



with the LHA. In the default case of a single collective coordinate only (K =1), the basic 
equations of this formalism can be written as 

V a = Xf 1 , (2.7) 

B*%ryfi = ^f!« ■ ( 2 - 8 ) 



(2.2) 

7i-=n 
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The second of these equation is the local harmonic (or local RPA) equation from which the 
method derives its name. The covariant derivative (denoted by ;) in the left-hand side of 
Eq. (ED is defined by 

v.* = v aP - r^y 7 , (2.9) 



with the affine connection T defined in the standard way in terms of the metric tensor B a 



■P 



as 



T/3 7 — 2~B a5 {B$p n + Bs 7) p — Bpy } g) . (2.10) 

The equations ( |2.7p and ( |2.8|) must be solved as a self-consistent pair, except at stationary 
points, where the RPA ( |2.8|) is independent of the force equation fl2.7|) . This allows us to 
bootstrap our way up from such a stationary point. This procedure constructs a path by 
finding successive points at which an eigenvector /' n of the covariant RPA equation (p.8[) 
satisfies the force condition ( |2.7| ) at the same time. It is worth noting that since the equations 
at every point can be solved independently, no computational error is accumulated while 
calculating the collective path, which is a problem with, e.g., the formalism of Goeke and 
Reinhard |14| . 



The quality of decoupling can be measured by comparing two different collective mass 
parameters that can be calculated in the theory. If we calculate the derivatives d^ a /dx l in 
terms of the tangents of the path, we find 

b - d ^u d ^ rain 

511 " d£ Baf> d£ ■ (2 - U) 

The other mass parameter can be obtained by using the eigenvectors f l a obtained from the 
covariant RPA equation, 

B n = fi a B aP fy (2.12) 

This is equal to (fin) -1 if the decoupling is exact. Therefore, we define the decoupling 
measure D as 

D = (-Bn) B 11 — 1 . (2.13) 

The size of this measure D indicates the quality of decoupling. The smaller its value, the 
better the decoupling. 



B. Removal of Spurious Modes 

In this section we discuss how to treat the spurious modes. A typical example is given by 
the Nambu-Goldstone (NG) mode associated with the violation of particle-number conser- 
vation. We have presented a method adding a constraint to the original LHA formalism, in 
order to find a collective subspace orthogonal to the NG modes ||. However, in the model 
to be discussed here this method does not work because of a divergence problem associated 
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with the fact that the model has an exact zero mass parameter, det(B al3 ) = 0. Instead, in 
this paper, we choose to remove the NG degrees of freedom explicitly. 

For the models to be discussed in the following sections, the modes associated with a 
change in average particle number are given by a linear combination of coordinates: 

n 

7 NG (0 = E c «r, (2.14) 

a=l 

where c a is a constant. The problem is that this mode leads to a zero eigenvalue of the 
mass, 

B af3 fJ' = B aP c p = . (2.15) 

This means that we cannot invert the mass matrix. The only sensible way to deal with this 
is to remove these degrees of freedom from our space, by defining a new set of coordinates, 
jtn _ These are required to satisfy 

B <*Pf» ^ o , for Va and y = 1, • • • , n - M , (2.16) 
B aP fy = , for Va and y = n - M + 1, • • - , n , (2.17) 

where we assume that there are M NG modes (y > n — M). Then, we may formulate the 
LHA in the space of n — M dimension, {^ M } M =i...., n -M, in which det(-B M ^) ^ 0. 

M V J% = (^)V,; , (2.18) 

m; = v:; = b^v vIii , (2.i9) 

where indices y, v, ■ ■ ■ are running only from 1 to n — M. Our aim is to provide a feasible 
method to calculate this LHA, namely to calculate the mass parameter B vv , potential V(£), 
and their derivatives. 

The second equation Q2.17 ) determines tangent vectors of the NG modes. The rest of 
coordinates / M for y = 1, • • • , n — M are arbitrary as long as their derivatives are linearly 
independent from the others. The full Jacobian matrix allows us to define the derivatives 
of inverse transformation, g a as the inverse of /, 

/5§J = #, f%~9l = V$- (2-20) 

Since all /£, are constant (independent of coordinates), all g 01 ^ are also constant and the 
derivatives f^p (or g^ a/3 ) all vanish. This implies that within the NG subspace the connection 
vanishes, T = 0, and the geometric character of the transformation of any tensor is fully 
determined in the subspace that does not contain the NG modes. One can use this to 
calculate the new mass parameter and its derivatives as 

= f%B<#f» p , (2.21) 
£y = ~f%B^fygl , (2.22) 

and the derivatives of potential as 

V, = ~g a »V a , (2.23) 
V,u = 9fflV a p . (2.24) 
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III. THE 0(4) MODEL 



We shall first study the properties of the single-shell 0(4) model. We define fermionic 
operators c] m and Cj m that create or annihilate a particle in the J z = m sub-state. In terms 
of these operators we define four pairing (P, P<, P, and P') and two multipole operators 
(N and Q) that close under commutation, and generate the algebra o(4), 



P^ c jm c jm J ^jm^jmPjm ' (^-l) 

m>0 m>0 

N = c im C im ) Q = &jmCj m Cj m , (3-2) 



J +1 for |m| < , , 

aim _ \ -1 for |m| > n/2 ' 

Here we need to require that the pair multiplicity Q = j + 1/2 is an even integer in order 
for the algebra to close. The sign of <7j m is chosen so as to mimic the behaviour of the 
matrix elements of the axial quadrupole operator (jm\r 2 Y 2 o\jm) , and we shall call Q the 
quadrupole operator in the remainder of this work, even though it does not carry the correct 
multipolarity. 

As is well-known, the algebra o(4) is isomorphic to su(2) © su(2). This can be made 
explicit in terms of the quasi-spin operators 



A+ = i(pt + pt) =Al , A = ^(N + Q-n) , (3.4) 
S + = |(pt-Pt) =B I, S = \ {N - Q - O) , (3.5) 
which generate two independent su{2) algebras, 



[A+,AJ\ =2A , \B + ,B_] = 2B , (3.6) 
[A , A±] = ±A± , [B , B ± ] = ±B ± , (3.7) 
[A^B^f] = . (3.8) 

The Hamiltonian of the model is chosen as a simple quadratic form in (some of) the 
generators of o(4), 

H = -GP ] P - ^kQ 2 , (3.9) 

and mimics the pairing-plus-quadrupole model that has been such a successful phenomeno- 
logical model in heavy nuclei ||. Even though the Hamiltonian looks simple, it does not 
have a closed- form solution (it does not have 0(4) dynamical symmetry). Nevertheless a 
numerically exact solution for the Hamiltonian ( |3.9| ) can be obtained by simple diagonalisa- 
tion. To this end one rewrites the Hamiltonian in terms of the quasi-spin operators A and 
B, 

H = —G (A + + B+) {A. + BJ) - 2k (A - B ) 2 . (3.10) 
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This Hamiltonian commutes with the total particle number N = 2(Aq + Bq) + Q, and there 
are no further constants of the motion. The pairing force tends to align the two quasi-spin 
vectors A and B, so as to obtain the maximal pairing matrix elements, while the quadrupole 
force tends to de-align them (to maximise (A Q — Bq) 2 ). In this picture, the non-integrability 
of the model, as well as the physics described, is related to the competition between the 
pairing and the quadrupole force. This is identical to a competition between alignment and 
de-alignment of the quasi-spins. 

For a fixed number of particles N = 2n , we construct from the vacuum state |0) all 
states with a constant number of generators A + and B + , 



\n , k a ) 



f -k a )\(%-n + k a )\ 



1/2 



Al a B 



no—ka 



10} 



(3.11) 



where < k a < n®. Finding the eigenvectors of the Hamiltonian now involves a trivial 
matrix diagonalisation in this basis of dimension (uq + 1). 



A. The Coherent-State Representation and the TDH(F)B Equations of Motion 



The mean-field description of the Hamiltonian ( |3.10| ) is most easily based on the use 
of a product of su{2) coherent states, one for the A^ sub-algebra, and another for the _B M 
one. Each of these states is characterised by a complex variable, z a and Zb ||15|| . The time- 
dependent mean-field dynamics in this parametrisation is the classical Hamiltonian problem 
we shall apply our methodology to. We can also parametrise the coherent state with four 
real angles P,|"i"5]j , 



\z a , z b ) = exp [z a A + - z* a A^ + z b B + - z^B_] |0) , 



' e x 

cos - cos — 
2 2 



fi/2 



exp 



e x 

tan - exp(— i<p) A + + tan — exp(— i?p)B-\ 
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(3.12) 
(3.13) 



where we have used 



2 



Zb 



-expHVO • 



(3.14) 



The time-dependent Hartree-Fock Bogoliubov (TDHFB) equations are in this case the 
classical equations of motion obtained from the stationary condition of the coherent-state 
action 5S = 0, where 



S[z] = J dt (z a , z b \id t - H\z a , z h ) 



e 



■ 2 X 



dt — I (J) sin — h ip sin — 



(3.15) 
(3.16) 



and 



T~L — (z a , z b\H\z a , z^ 



(3.17) 



S 



In order to facilitate our work we introduce real canonical variables £ Q and ir 



n 



n 



n 



n 



£ = — sin \z a \ = — sin - , £ = — sin 1 2^ | = — sm — 



(3.18) 

7Ti = arg^ ; = —<t> , n 2 = aig{z b ) = -ip . (3.19) 
Since these variables are canonical, the equations of motion are of Hamiltonian form 



2 2 ' * 2 
7T 2 = arg(z 6 ) = -ip 



_&h , a _m 



(3.20) 



where the classical Hamiltonian ( |3.17|) is the coherent state expectation value of the Hamil- 
tonian rewritten in terms of canonical variables, 



Ti = Tip + T~L 

G 
16 



+ 32n- 1 X)« 



»\2 



(3.21) 
(3.22) 

(3.23) 
(3.24) 



The terms proportional to Q^ 1 in Eqs. (|3.22|) and ( |3.23| ) originate from the exchange (Fock) 
contributions. 

Expanding the Hamiltonian with respect to n up to second order, the mass parameter 
B a/3 (£) and the potential V(£), which both depend only on the coordinates, are defined as 



a/3 



— 5 a p ( S a ^2 S~A — S a Sj3 , (no summation with respect to a) 
8 L V 7 / 
v(0 = V P (0 + V Q (£) , 

^(o = -2*j(E(-r +1 ^) +^ l En^-nj • 



(3.25) 

(3.26) 
(3.27) 
(3.28) 

(3.29) 



In the rest of this work we shall invoke the Hartree-Bogoliubov approximation, and we 
neglect the exchange terms (Fock terms) in Eqs. (|3.28| ) and (|3.29|) . 



B. Requantisation 

In this section, we discuss the problem of defining a requantisation procedure and the 
consequences of the adiabatic truncation with respect to momentum. The classical limit of 
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the single- j Hamiltonian has two constants of motion: 7i = E and (N) = 2 J2 a = -Wo- 
Since the phase space is four dimensional, this implies the complete integrability of the 
system, and there is a two dimensional torus on which all classical orbits lie. Due to this 
special feature of this model, one can apply the Einstein-Brillouin- Keller (EBK) quantisation 
condition. This has been done in Ref. |J and good agreement with the exact results has 
been obtained for both energy spectra and transition amplitudes. However, it is impossible 
to extend this quantisation method to non-integrable systems like the ones we will discuss 
in the following sections. We wish to use the same quantisation procedure for the simplest 
form of the model and the more complicated cases discussed later on, and shall turn to our 
favourite technique first. 

After truncation of the Hamiltonian up to second order in momentum, we can define 
a collective Hamiltonian by evaluating its value for points on the collective space £ which 
is parametrised by x 1 and p 1 (strictly this is the co-tangent bundle over £), since we have 
chosen x a = p a = 0; a = 1, • • • , n — 1, 

n col = K d | E = ^B l \x l )pl + V{x l ) , (3.30) 



^V) = E/i/i^ Q/3 (r = g a (x\x a = o)) , (3.31) 

a/3 

Vi^x 1 ) = v(i a = g a (x\ x a = 0)) . (3.32) 



Since the scale of collective coordinate x l is arbitrary, we choose to normalise f \ so as to 
make B 11 = 1. Subsequently, the Hamiltonian H co i is quantised in this flat space as [Q] 

Hcoi = -\^ + V(x), (3.33) 

where we employ the simplification of denoting (x l ,pi) by (x,p). 

In order to evaluate the matrix elements of a one-body operator F (either diagonal 
or transition matrix elements), we first obtain the collective classical representation of the 
operator F, which in keeping with the adiabatic approximation is expanded in powers of 
momentum, 

oo oo 

F(x,p) = ^,7r)| s = (z|F|s)| E = £ ^(e,vr)| s = £ • (3.34) 

i=0 i=0 

Here J 7 ® is the term of i-th order in it. The function T is requantised, by making the 
replacement T{x,p) — > F(x, 4M, at which point one will have to confront the problem of 
operator ordering between x and p. We shall avoid this problem by keeping, invoking once 
again the assumption of slow collective motion, only the zeroth order term j 7 ^. It is an 
interesting question what the effect of higher order terms will be. This is clearly outside the 
scope of the present work, and requires further investigation. Fortunately, in the current 
model, we have no ambiguity for the quadrupole operator Q because QW = o for i ^ 0. For 
convenience we denote the classical limit of the quadrupole operator by q. The transition 
matrix elements can thus be calculated by the one-dimensional integral, 

/ dx^* n ,(x)F ( -°\x)^ n (x) , (3.35) 
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where \l/ n are the eigenf unctions of the collective Hamiltonian fl3.33|) , 

From the number of coordinates and momenta found (2+2) we see that the configuration 
space of the single- j shell model is two-dimensional. Since there is a zero mode (det(B a @) = 
0) corresponding to the NG mode associated with the particle number violation, one may 
obtain a one-dimensional path S by simply applying the prescription in Sec. [11 B| (without 
the application of LHA)! Rather than plotting this path we have chosen to represent the 
results of requantisation for energies and transition strengths. These results are presented 
in Fig. |l|. We obtain good agreement with the exact results over a wide range of parameters 
except very close to a pure quadrupole force. Due to the peculiar nature of the quadrupole 
operator the mass parameter goes to zero in this limit (G = 0), and there is no kinetic 
term. Thus an eigenstate of Hamiltonian is a coordinate eigenstate \x) at the same time. 
Then, the periodic nature of the momentum becomes important, which we have ignored in 
our calculations. Taking account of the periodicity of momentum, one finds that coordinate 
operator x should have discrete eigenvalues. In order to check that it is possible to deal 
with this problem, we have expanded the Hamiltonian up to second order with respect to 
the coordinates rather than momenta, keeping all order in momenta. We have also imposed 
periodic boundary conditions on the wave function ^f(p) = \J/(p + 7r/2). The result of this 
quantisation is shown in Fig. ^. The agreement is good in the no-pairing limit G — 0, 
while it is not as good as the standard quantisation anywhere else. Since we are not really 
interested in the (integrable) pure quadrupole model, but rather in competition between the 
pairing and quadrupole forces, we shall ignore the G = limit in the rest of this work. We 
shall thus follow the conventional adiabatic quantisation procedure in coordinate space as 
described above. 



IV. THE MULTI-0(4) MODEL 

A. The Model and the Hamilton's Form of TDH(F)B Equation 

It is a greater challenge to our approach to study the multi-shell case. There exists a 
straightforward extension of the model, by addition of the individual pairing generators, and 
summing the quadrupole operators of each shell with a weight factor (we shall thus not have 
a direct coupling between the different shells). The operators are, for j = jx,j2, • • • , Ja, (for 
each shell ji we take the pair degeneracy f2j = ji + 1/2 to be even) 

= E C }im c ] t fh t ) = E a jimi c j i m i c jifhi ) (4-1) 

i,m,i>0 i,m,i>0 

N — '^^ C ] i m i C jim i j Q = '^^Qi°'jimi C } i m i C jimi (4-2) 

irrii irrii 

where represents the magnitude of quadrupole moment carried by single-particle states. 
For each shell we can define quasi-spin [su(2) © su(2)] generators A* and B* in a manner 
similar to Eqs. fl3.4| ) and ( ^.5| ). We choose a slightly more general Hamiltonian than in the 
previous chapter by adding a term containing spherical single-particle energies, 

H = E tjcLcjm - GP^P - ]-kQ 2 . (4.3) 

jm 
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The exact solution can be also obtained by diagonalisation in a basis set 



®i*s°,*? ) > = nw fci ( s ;)* i i o >. 



(4.4) 



i=i 



i=l 



where < k l a , k\ < Q{/2 and ^Zi{k l a + k l b ) = tiq = Nq/2. This is no longer as trivial a 
calculation as before, since the dimension of the basis increases rapidly with the number of 
shells A, but can still be done, provided that one chooses A sufficiently small. On the other 
hand, since the dimension of TDH(F)B configuration space increases only linearly with A, 
the amount of effort required for the ALACM calculation is still rather small. The time- 
dependent mean-field state is obtained through the use of the coherent-state representation 
as before, and is given by 



exp 



E4^ 

Li=l 



z?AL + z\B\ - z?BL 



|0>, 



=n 



Xi 



cos — cos — exp 
2 2 



tan exp(-i0i)A l + + tan y exp(-iipi)B l + 



(4.5) 
(4.6) 



A possible choice of canonical variables is 

for a — i — 1, • • • , A , 



r 



Tin 



T sin 2 f , 

*f sin 2 f , for a = A + i = A + 1, • • • , 2A , 

— 0i , for a = i — 1, • • • , A , 

-ipi , for a = A + i = A + 1, • • • , 2A . 



(4.7) 
(4.8) 



It is convenient to allow the indices of e, q and Q to range from lto 2A by copying the 
original list of parameters, e.g., 



J e,i (i — a) for a — 1, • • • , A , 
6q ~ \ a (i = a - A) for a = A + 1, • • • , 2A . 

Using these definitions, the classical Hamiltonian can be given in the compact form 

7~t = h sp + Hp + TLq , 



G_ 
16 



+32£^(n 2 



0"„ = 



2^(0* - 2^) , 

f +1 for a = 1,---,A , 
1-1 for a = A + 2A . 



(4.9) 

(4.10) 
(4.11) 

(4.12) 

(4.13) 
(4.14) 



The adiabatic limit of this Hamiltonian is 
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ai3 



B a(3 



G 



a/3 



S a S-y ) S a Sfi 



(4.15) 

(4.16) 

(4.17) 
(4.18) 



Once again, the terms in Eqs. ( |4.11| ) and ( (4.12[ ) proportional to Q a 1 arise from the exchange 
contributions, and will be neglected. 



B. Results for Transitional Nuclei 

Using the LHA we identify the collective degree of freedom amongst the 2A coordinates. 
As before we first must remove the NG mode corresponding to a change in particle number 
explicitly, due to the zero mass parameter associated with this mode. The particle number 
is simply given by the sum of the numbers for the individual shells, M = 2 J2 a £ a - ft is easy 
to show that 

]T B^Np = 2 Ba(3 = • ( 4 - 19 ) 

(3 (3 

Thus we apply the prescription discussed in Sec. [II B| to this model, and use the LHA 
to determine the collective path in the remaining (2A — l)-dimensional coordinate space. 
Since we shall mainly investigate how the LHA can deal with the transition spherical to 
deformed, it is sufficient to study only two shells. We take the size of these shells to be 
equal, Qi = Q 2 = 10, and put 16 particles in the available space. We split the degeneracy 
by taking e± — and e 2 = 1, and we use a different value of the "quadrupole moment" for 
each shell as well, q± = 3 and q% = 1. The pairing force is fixed at G = 0.3, and we only 
vary the quadrupole force strength k. 

We show a representation of the collective potential energy and the collective coordi- 
nate in Fig. [| We represent these quantities as a function of the expectation value of the 
quadrupole operator. As an alternative to the LHA machinery, we have also performed a 
simple constrained Hartree-Bogoliubov (CHB) calculation, where one imposes a value for 
the expectation value of the quadrupole operator. We determine the mass for this case by 
replacing the RPA eigenvector by the coordinate derivative of the quadrupole expectation 
value. We then renormalise the coordinate to obtain a constant mass. 

We have investigated a full shape transition scenario, where we have changed the strength 
of the quadrupole interaction so that the collective Hamiltonian changes from spherical and 
harmonic for case (a) to flat for case (b) to deformed for case (c) . We see that the difference 
between the LHA and CHB calculations is relatively small. This is also borne out by 
the spectra and transition strength in Fig. [|. We can see the similarity between the two 
approximate calculations. If anything, the LHA gives slightly better results than the CHB 
based calculations. We seem to be unable to reproduce the large density of states found in 
the exact calculation for "deformed" nuclei, where there are indications from the transition 
strengths that some of the states in the approximate calculations are split into several of 
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the exact states. Note, however, that at an excitation energy of 6, we are 5 units above 
the barrier, so this may just be due to our choice of parameters. The shape mixing in the 
low-lying excited states appears to be described sensibly, however. We would have liked 
to be able to choose an even large value of k, but if we do that the system collapses to 
the largest possible quadrupole moment in the model space, which leads to all kinds of 
unphysical complications. 



V. A MULTI-SHELL 0(4) MODEL WITH NEUTRONS AND PROTONS 

In heavy nuclei, the number of neutrons is normally (much) larger than that of protons, 
which often leads to a radically different shell structure at the Fermi surface for neutrons and 
protons. In order to perform a model study of such phenomena, where we can still perform an 
exact calculation, we adapt the multi-shell 0(4) model introduced in the previous section to 
one describing systems with both neutrons and protons |7[. We shall then use this model to 
analyse the collective dynamics of shape-coexistence nuclei, as observed for instance in semi- 
magic nuclei. At the same time we shall concentrate on the diabatic/adiabatic dichotomy 
already mentioned in the introduction. 

The model is a simple extension of the multi-shell 0(4) model in the previous section, 
with the main difference that we do not have pairing between proton and neutron orbitals, 

H = H n + H p + H np , (5.1) 

H n = ^ €iC^ jmj Cj imi — G n PlP ri — ^Qli j (5-2) 

i£n,m,i 

Hp = ^i C jim,i C jimi ~ GpPpPp — —K-Qp , (5-3) 

H n p = nQnQp i (5-4) 

where P n ( P ) and Q n ( P ) are the pairing and quadrupole operators for neutrons (protons) (see 
the definitions in Sec. [IVD . In this model there are two trivial NG modes associated with 



the change of neutron and of proton number, which can both be removed explicitly in the 
manner discussed before. 

We study this model for a single-shell for neutrons, with pair degeneracy Q n = 50, 
containing 40 particles. We take the single-particle quadrupole matrix element q n = 1, 
and use a pairing strength G n = 0.3, and assume zero single-particle energy. For protons 
we take two shells, both with Q pl = Q p2 = 2, q p i,q P 2 = 2, having single-particle energies 
e p i = —e p2 = 5. We study two different set of interaction parameters, both with k = 0.1. 
The first is G n = G p = 0.3, the second has the same neutron pairing strength, but G p = 10. 

The collective potential energy for the weaker pairing strength shows a very interesting 
behaviour, with two shoulders appearing in the CHB collective potential energy. This is 
what is normally called the adiabatic potential energy, and the shoulders arise from an 
avoided crossing. As in our previous work || the adiabatic LHA method chooses a diabatic 
(crossing) set of potential energy curves. These shape-coexistence minima are related to 
2p-2h excitation in the proton model space, promoting two particles from the lowest Nilsson 
orbitals to the down-sloping ones. This is of course very similar to the phenomena observed 
in shape coexistence in semi-magic nuclei. 
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We get another surprise for the strong pairing case. Here the collective potentials look 
very similar, and rather structureless, but the collective coordinates are different. This can 
be traced back to the fact that the collective coordinate in the LHA is not q p + q n , but a 
different combination. At the minimum, the lowest RPA mode correspond approximately 
to q n + j$q P - This is similar to the situation analysed in great detail in Ref. fllEfl , and once 
again shows the importance of self-consistency in the selection of the collective coordinate. 
The real collective coordinate is not the mass-quadrupole operator! 

In Fig. ^| we show the consequences of these results for the requantisation. For the weak 
pairing case the diabatic picture obtained through the LHA gives almost perfect results, 
whereas the CHB potential energy curve fails to provide the correct answer. In the case of 
the strong pairing the incorrect choice of the collective coordinate leads to too large a level 
spacing in the CHB calculations, whereas the LHA and exact calculations agree again. 

For the weak pairing case the exact calculation almost exhibits the doublet structure 
found from the diabatic potential energy curves; the splitting is less than one part in 10 5 
in the exact calculation. Nevertheless, the exact calculation consists of symmetric and 
antisymmetric states, which leads to the transition matrix element 22.3 between these two 
states. Since in the LHA calculation the states do not mix, we have printed the diagonal 
matrix element instead. For very weak mixing this is the correct comparison, as is borne 
out by the results. 

The decoupling measure D, Eq. ( |2.13| ), is small for all these cases. The worst case is the 
strong pairing case, where it rises from for q = to 3 x 10~ 3 for q = 40. For this reason 
we also believe that the "scalar Berry potential" p| will be small, and we have not included 
this, or any other quantum corrections, in our calculations. 



VI. CONCLUSION 

In the area of shape coexistence, as studied by the model calculations discussed in this 
paper, it seems that on the whole the LHA is an excellent method to obtain sensible results, 
with reasonable effort. Even though the size of the space is much smaller than that used in a 
straightforward diagonalisation, the use of the local RPA, which needs to be performed many 
times at every point along the collective path, may seem prohibitive for realistic calculations. 
Fortunately, many problems can be studied with simplified separable forces, as in the model 
discussed here. We are at the moment considering the old pairing-plus-quadrupole model, 
that has been applied so successfully in the physics of heavy nuclei. Such a model can 
be dealt with much more straightforwardly than more microscopic Skyrme or Gogny-force 
based approaches. This will allow us to shed light on a good treatment of shape-coexistence, 
as well on the interesting question on the choice of the collective (cranking) operator, which 
was already found to be non-trivial in certain limits of the 0(4) model. 

One might argue that even that is not enough, and we should really adopt the full horror 
of standard nuclear many-body approaches. We believe that we can address this problem, 
and are actively considering the approaches available to us. As must be obvious from the 
discussion given above, an efficient calculation hinges on efficient solution of the RPA. We 
are investigating two approaches to this problem, the use of iterative diagonalisation of the 
RPA using Lanczos procedures, or the approximation of the RPA by using separable forces 
|17| |, which can be diagonalised much more efficiently. 
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In the present work we have not included any quantum corrections, nor the "scalar Berry 
potential" arising from choosing a slow degree of freedom. In our previous work || we have 
shown that these corrections tend to improve the spectrum, so that the results presented 
can even be improved slightly. We are strongly encouraged by the ability of our method 
to choose the right coordinate and the correct potential energy curves for problems that 
involve pairing and "quadrupole" degrees of freedom, and we expect results for realistic 
nuclear models in the near future. 

ACKNOWLEDGMENTS 

This work was supported by a research grant (GR/L22331) from the Engineering and 
Physical Sciences Research Council (EPSRC) of Great Britain. 



16 



REFERENCES 



[1] A. Klein, N.R. Walet and G. Do Dang, Ann. Phys. 208, 90 (1991). 

[2] T. Nakatsukasa and N.R. Walet, Phys. Rev. C57, 1192 (1998). 

[3] R. Piepenbring, B. Silvestre-Brac and Z. Szymanski, Nucl. Phys. A348, 77 (1980). 

[4] K. Matsuyanagi, Prog. Theor. Phys. 67, 1441 (1982); Proceedings of the Nuclear Physics 
Workshop, Trieste, 5-30 Oct. 1981. ed. C.H. Dasso, R.A. Broglia and A. Winther (North- 
Holland, 1982), p.29. 

[5] M. Matsuo, Prog. Theor. Phys. 76, 372 (1986). 

[6] T. Suzuki and Y. Mizobuchi, Prog. Theor. Phys. 79, 480 (1988); Y. Mizobuchi, Prog. 

Theor. Phys. 65, 1450 (1981). 
[7] T. Fukui, M. Matsuo and K. Matsuyanagi, Prog. Theor. Phys. 85, 281 (1991). 
[8] P.O. Arve and G.F. Bertsch, Phys. Lett. B215, 1 (1988). 
[9] D.R. Bes and R.A. Sorensen, Adv. Nucl. Phys. 2 (1969) 129. 
[10] T. Kishimoto and T. Tamura, Nucl. Phys. A192, 246, (1972); Nucl. Phys. A270, 317, 
(1976). 

[11] J.L. Wood, K. Heyde, W. Nazarewicz, M. Huyse and P. van Duppen, Phys. Rep. 215, 
101 (1992). 

[12] R. Bengtsson and W. Nazarewicz, Z. Phys. A334, 269 (1989). 
[13] D.L. Hill and J.A. Wheeler, Phys. Rev. 89, 1102 (1953). 
[14] PG. Reinhard and K. Goeke, Rep. Prog. Phys. 50, 1 (1987). 

[15] A. Perelomov, Generalized Coherent States and Their Applications, (Springer- Verlag, 
Berlin, 1986). 

[16] G. Do Dang, N.R. Walet, and A. Klein, Phys. Lett. B322, 11 (1994). 
[17] V.O. Nesterenko, W. Kleinig, V.V. Gudkov, and J. Kvasil, Phys. Rev. C 53, 1632 
(1996). 



17 



400 



300 




(2) adiabatic-2 



400 ■ 



300 ■ 



200 • 



100 




400 



300 



200 



100 



1 (3) exact 



\l 



,19.9, 
,19.3 

,16.8 



(a) (b) (c) 



(a) (b) (c) 



FIG. 1. The excitation energies and transition matrix elements |(n — l|Q|rc)| in the single- j 
shell 0(4) model as a function of the ratio between the strength of the pairing and that of the 
quadrupole force. The case (a) is a weak quadrupole force, 2k/ G = 0.079, (b) a medium sized one, 
2k/ G = 1.63 and (c) a very strong one, 2k/ G = 12.7. The three panels give our standard adiabatic 
quantisation (1), the results from the "adiabatic in coordinates" method (2), and the exact results 
(3). 
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Tt/2 

tan" 1 (2K/G) 

FIG. 2. The excitation energies in the single- j shell 0(4) model for 40 particles in a shell with 
0, = 100 as a function of the parameter essentially the ratio between the strength of the pairing 
and that of the quadrupole force. The left end corresponds to the case of pure pairing force and 
the right end to the pure quadrupole force. The solid lines are the exact results, and the dotted 
lines represent the standard requantisation of the adiabatic collective Hamiltonian. The dashed 
lines represent the expansion in terms of coordinates discussed in the text. 
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FIG. 3. The collective potential energy V(q) and the collective coordinate x (normalised to 
unit mass) as a function of the quadrupole moment q = (Q). We show both the LHA (dashed 
line) and the CHB (solid line) results in each figure. These results are, for this model, essentially 
indistinguishable. The case (a) corresponds to a weak quadrupole force (k = 0.01), (b) to a slightly 
stronger force (k = 0.03) and (c) to the strongest, k = 0.35. 
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FIG. 4. The excitation energies E x and transition matrix elements |(n'|Q|n)| (numbers next to 
arrows) in the two-shell case discussed in the text. The three cases (a), (b) and (c) correspond to 
those discussed in Fig. ^. In each of the three panels the left one (1) is obtained from requantising 
the CHB, the middle one (2) is obtained from requantising the LHA result, and the right one 
corresponds to exact diagonalisation. The dashed line shows the lowest RPA eigenvalue. 
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FIG. 5. The collective potential energy V(q) and the collective coordinate x as a function of 
the quadrupole moment q = q p + q n . We show both the LHA (dashed lines) and the CHB results 
(solid line) in each figure. The case (a) corresponds to a weak proton pairing force (G p = 0.3), (b) 
to a strong pairing force (G p = 10). The rest of the parameters are given in the main text. 
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FIG. 6. The excitation energies E x and transition matrix elements |(n'|Q|n}| (numbers next to 
arrows) in the proton-neutron case discussed in the text. The two cases (a) and (b) correspond 
to those discussed in Fig. ||. In each of the panels the left one (1) is obtained from requantising 
the CHB, the middle one (2) is obtained from requantising the LHA result, and the right one (3) 
corresponds to the exact diagonalisation. Levels denoted by thick lines are doubly degenerate. The 
numbers next to the arrows denote teh size of the transition matrix elements. 
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